Vacuum local and global electromagnetic self-energies for a point-like and an extended 
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We consider vacuum electric and magnetic energy densities (or equivalently field fluctuations) in 
the space around a point-like or an extended ground-state field source, and discuss the problem of 
their singular behavior at the source's position. We show that the assumption of point-like source 
leads to a divergence of the renormalized energy densities at the position of the source and analyze 
in detail the structure of such singularity in terms of a delta function and its derivatives. We then 
consider the case of an extended source, smeared out over a finite volume and described by an 
appropriate form factor, and show that in this case all divergences in the local energy densities are 
removed. Our results for the structure of the divergences in the field energy densities also resolve 
an apparent inconsistency between the vacuum expectation value of the field Hamiltonian and of 
the field energy densities integrated over all space. 

PACS numbers: 12.20.Ds, 03.70.+k 



I. INTRODUCTION 

Vacuum fluctuations and the existence of the zero- 
point energy of the electromagnetic field are a remarkable 
prediction of quantum electrodynamics [l|. They have 
been extensively investigated in the literature, especially 
in relation with Casimir and Casimir-Polder forces, which 
are long-range electromagnetic interactions between neu- 
tral macroscopic objects (metals or dielectrics), between 
atoms and surfaces or among neutral atoms or molecules 
0-01 ■ The Casimir effect, in particular, is related to the 
modification of the zero-point energy as a consequence of 
a change of the boundary conditions on the electromag- 
netic field 0. 

Despite of its oddity, the idea of a force generated by a 
change of vacuum fluctuations in the presence of bound- 
ary conditions has stimulated great interest in the liter- 
ature, both from theoretical and experimental point of 
view [||. These purely quantum effects have now been 
experimentally confirmed with remarkable accuracy @- 
8]. However, there are many aspects concerning with 
vacuum fluctuations and vacuum energy that are still 
unclear. 

A controversial issue concerns with the appearance of 
surface divergences (and their cut-off dependence) in the 
calculation of field energy densities, in the presence of 
metallic boundary conditions [9l4l3j . The physical origin 
of these divergences has been recently questioned in the 
literature and the possibility of removing them through 
a suitable regularization procedure has been discussed in 
the case of a scalar field [ll|, [l3[ • Generally speaking, the 
appearance of surface divergences in the field energy den- 
sity is not surprising: boundary conditions are in fact a 
convenient oversimplification of the interaction between 
matter and fluctuating fields. An ideal boundary con- 
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straints all field modes, at any wavelength and this gives 
rise to singular energy densities on the boundary. For 
example, it is well-known that, in the case of a perfectly 
reflecting plate, renormalized electric and magnetic fluc- 
tuations or energy densities (i.e. after subtraction of the 
homogeneous energy density existing also in absence of 
the boundary) (E 2 ) and (B 2 ), diverge at the vacuum- 
conductor interface. These divergences can be removed 
by introducing an appropriate exp onential upper cut-off 
in the frequency integrations [l2l FLlj . It has been also 
shown that they disappear in the case of a boundary with 
a fluctuating position [ll|, and it has beenguessed that 
imperfect boundaries could remove them [la ]. 

Recently, the electromagnetic field fluctuations near a 
dielectric-vacuum interface have been investigated [HI, 
[HI ] . The structure of the surface divergences in the limit 
of a perfect conductor has been discussed in detail [ilj . 
Also, local energy densities and surface divergences have 
been explored near fluctuating boundaries [llj or in the 
vicinity of a soft wall modeled by a potential which is 
monomial in z, z being the distance from the wall [17] . 

In spite of the increasing interest on the subject, many 
topics are not yet well understood. For example, the pres- 
ence of divergences in the field energy density presents 
serious problems when the coupling to gravity is consid- 
ered, because considerable gravitational effects should be 
observed [13, Ha]. In fact, even if considering real surfaces 
might eliminate the divergences associated to idealized 
boundaries, energy densities at surfaces and their effects 
could be large and should be carefully considered |l9Ll20j. 
Also, it should be stressed that electromagnetic energy 
densities can be experimentally investigated through the 
Casimir-Polder interaction energy with appropriate elec- 
trically and magnetically polarizable bodies [21| . 

All the considerations above make relevant investigat- 
ing the structure and the physical origin of divergences 
of the vacuum energy densities in the presence of bound- 
aries or polarizable bodies. In [l4| we have investigated 
the origin and structure of surface divergence at the in- 
terface between an ideal conductor and the vacuum space 
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exploiting an appropriate limiting procedure from a di- 
electric to a metal; this has also allowed us to show that 
surface divergences in the electromagnetic energy density 
are essential in order to have consistency between inte- 
grated local energy densities and global field self-energies. 

It is thus natural investigating the case of renormalized 
vacuum electromagnetic energies near an atom or a po- 
larizable body, and this will be the main subject of this 
paper. 

In this paper we study the electromagnetic energy den- 
sities (or, equivalently, the held fluctuations) surrounding 
a dressed point-like and a dressed extended source of the 
electromagnetic held, whose properties are given by its 
dynamical polarizability and its form factor; in partic- 
ular, we shall concentrate on their singular behavior at 
the source's position. Compared to previous works on the 
subject mentioned above, in the present case we have a 
dynamical source of the electromagnetic held in place of a 
boundary condition. Some aspects of the vacuum fluctu- 
ations of the electromagnetic held in presence of a point- 
like source have been already explored in the literature 
(see for example [22[ and references therein), as well as 
its connection with Casimir-Polder interactions between 
atoms and/or macroscopic bodies (23j . but its singular 
behavior was not considered. It is known that, under ap- 
propriate conditions (in the so-called far zone), both the 
renormalized electric and magnetic energy densities sur- 
rounding a point-like ground-state atom are proportional 
to 1/r 7 , r being the distance from the atom. They clearly 
diverge at r = 0, as well as their sum and their integral 
over all space. On the other hand, the renormalized total 
held energy (expectation value of the held Hamiltonian, 
after subtraction of the zero-point terms) can be shown 
to be zero (see Sec. [TTJ) . This shows an inconsistency be- 
tween local and global self-energies, analogous to the case 
of the ideal metallic boundary condition discussed in [l4[ , 
where it is shown that only a careful consideration of sur- 
face divergences at the metal-vacuum interface permits 
to restore a full consistency between the two approaches. 
This gives a strong indication that also in the case of a 
point-like atom there should be extra singular terms at 
the atomic position. It has been argued that the singu- 
lar behavior arises from the oversimplified assumption of 
point-like source and dipole approximation 2J| . In some 
sense, point-like sources are conceptually similar to ideal 
boundaries, so singular energy densities at the position 
of the source are to be expected. 

These considerations suggest that divergences could be 
removed with a suitable regularization procedure involv- 
ing the use of an appropriate upper-frequency cut-off in 
the calculations. In this paper we shall explore an al- 
ternative procedure to treat singularities in the energy 
densities at the position of a field source, showing the 
presence of extra singular terms localized at r = 0, in 
addition to the 1/r 7 term. This will allow us to prove 
that there are not inconsistencies between global and lo- 
cal held self-energies, provided the energy-density diver- 
gences localized at the position of the source are prop- 



erly included. We shall also consider an extended source 
of hnitc dimension interacting with electromagnetic field 
in the vacuum state. The source is modeled as a col- 
lection of elementary neutral sources, smeared out over 
a hnitc volume with a density described by a function 
p(r), whose Fourier transform gives a form factor for the 
source. This form factor, which appears in the matter- 
held interaction Hamiltonian, plays the role of a regular- 
ization factor in the source-field interaction, and it pro- 
vides a finite length-scale for the source-field coupling 
which cuts off the contribution of high frequency modes. 
This ensures that the local energy densities of the elec- 
tromagnetic held are well-defined and finite everywhere, 
thus eliminating the problem of divergences present in 
the case of the point-like source. Global and local self- 
energies are fully consistent in this case, even for a van- 
ishing size of the source. The case of a point-like source 
can be then obtained by an appropriate limit procedure. 
This result also indicates that the procedure adopted is 
equivalent to consider an ideal point-like source together 
with a suitable physical cutoff as a model to describe 
realistic situations [19J. 

The paper is organized as follows. In Sec. [II] we con- 
sider a point-like source in its dressed ground-state in- 
teracting with the quantum electromagnetic field in the 
vacuum state. We assume the source located at r = and 
investigate the local and global (i.e. integrated over all 
space) properties of the electric and magnetic energy den- 
sities in the vacuum space surrounding the source. We 
show that the assumption of a point-like source leads to 
divergences in the local electromagnetic energy densities, 
and obtain their explicit mathematical expression. We 
show that a singular behavior at r = is at the origin 
of the apparent discrepancy between the total electro- 
magnetic self-energy obtained as the expectation value 
of the held Hamiltonian and as a space integral of the 
held energy density. We hnd that such inconsistence is 
completely removed if the singularity at the source posi- 
tion r = is correctly taken into account, by virtue of a 
subtle cancelation between diverging quantities at r = 0. 
In Sec IIIII we consider the case of an extended source with 
a finite size and analyze in detail the behavior of (local 
and global) field energy densities in the space around 
the source. We discuss how the proposed model of ex- 
tended source eliminates divergences and singularities in 
the electric and magnetic energy densities, and also solve 
the mentioned inconsistencies in the calculation of the 
total electromagnetic energy for any nonvanishing size of 
the source. Sec. HVl is finally devoted to some concluding 
remarks. 



II. ELECTRIC AND MAGNETIC ENERGY 
DENSITY AROUND A POINT-LIKE SOURCE 

In order to discuss the problem of the divergences of 
the (renormalized) electromagnetic energy densities, we 
first consider the case of a point-like atom in its ground- 
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state interacting with the electromagnetic field in the vac- 
uum state. We suppose the atom located at r = and 
we describe the atom-field system adopting an effective 
Hamiltonian obtained from the multipolar form of the 
atom-field Hamiltonian in dipolc approximation through 
an appropriate transformation pH | 



H = H A + H F - a ( fc ) d kj(0) • d(0) 



(1) 



kj 



(we are working in the Coulomb gauge) . Ha is the atomic 
Hamiltonian and Hf is the electromagnetic field Hamil- 
tonian; d(r) is the transverse displacement field, whose 
Fourier components are dkj(r): outside the atom, it co- 
incides with the total (transverse plus longitudinal) elec- 
tric field. Introducing the usual annihilation and creation 
bosonic operators auj and a fc ■, the Hamiltonian of the ra- 
diation field assumes the well-known expression 



HF = hj (d2(r) + b2 ( r )) d3r = E h ^ a h a ^ 

kj 



(2) 



(in the last equality the zero-point energy has been sub- 
tracted), with 

d(r) =^Y: (^) 1/2 (^^^' r - e^t.e-*--) (3) 



kj 



b(r) =*£ (^) 1/2 (S ki a kj e^ r -S^.at . e -*-)(4) 



kj 



where e^j are the polarization vectors and b^j = e^j x 
k. Finally, the atomic isotropic dynamical polarizability 
a(k) is given by 



2 x 

a(k) = - 



1 

mg 
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(5) 



We now evaluate the expectation values of the squared 
electric and magnetic fields, on the dressed ground-state 
©. After straightforward algebraic calculations we get, 
at order a(k), 
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(0|W rf (r) | 0) = -^(0 | d 2 (r) | 0) 

kk' 



k + k' 



kj k'j' 

x (a(k) +a(fc'))e i(k+k,) ' r + c.c., 
(0|W m (r)|0) = ^(0|b 2 (r)|0) 

= ^ (7) EEfe • g k'j')(^kj ' bk'j') 
kj k'j' 

x (a(fc) +a(fc'))e i(k+k,) ' r + c.c., 



(7) 



kk' 



k + k' 



(8) 



where we have subtracted the spatially uniform zero- 
point contributions, thus obtaining the renormalized en- 
ergy densities (from now onwards, we will consider only 
renormalized self-energies). The procedure adopted so 
far is rigorous up to terms of order a(k). 

In order to calculate explicitly the quantities and 
([8]) , we approximate the dynamical polarizability with its 
static value. This is equivalent to consider the so-called 
far-zone (r > c/uiq), where ujq is a main characteristic 
atomic frequency: in this zone the prevailing contribu- 
tion to the energy densities is given by virtual photons 
with frequencies u>k < ^0 [HJ . Thus we can approximate 
ujq + ujk ~ luq, so that the dynamic polarizability a(k) 
reduces to the static one, a — a(0). This approximation 
is mathematically consistent for any distance r, provided 
ujq is sufficiently large, and it will make our considera- 
tions independent from a specific structure of the field 
source. 

In the continuum limit (V — > 00), after sums over po- 
larizations and angular integrations, we get 



where fx mg is the matrix element of the dipole moment 
operator between the atomic state m and the ground 
state g, and E mg — E m — E g . We are interested in inves- 
tigating local and global properties of the energy density 
of the virtual electromagnetic field surrounding the point- 
like source in its dressed ground state. In particular, we 
shall calculate the quantum average of the electric and 
magnetic energy densities on the dressed ground state, 
that is the lowest-energy eigenstate of the total Hamilto- 
nian (fl}. Using Rayleigh-Schrodinger perturbation the- 
ory, the dressed ground-state is easily obtained at first 
order in a(k) 
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' ~vl^l^ a(fc) k + k' 6kj ' 6k ' J 



V 

9, Ikjlk'j')) 



kj k'j' 



(6) 



where g denotes the atomic ground state. 



(0\H e i(r) I °> = 4^3 / dk I dk ' 



00 />oo 



xQ E {k,k' ,r] 



k 3 k' 3 
k + k' 



he 



(0 I H m (r) I 0) = - — / dk I dk' 



OO /'OO 



xQ M {k,k\r 



k 3 k' 3 
k + k r 



(9) 



(10) 



where QE(k,k',r) and Qu{k,k',r) are expressed in 
terms of spherical Bessel functions as 

QE(k,k',r)=j (kr)j (k'r)-j (kr) J -^- ) 
J]il ' r) 3o(k , r) + -^ J1 (kr) J1 (k'r) (11) 



kr kk'r 
Q M (k,k',r) = 2j 1 (kr)j 1 (k'r) 



(12) 
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Using the relation 
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k + k' 



-(k+k') v 



(13) 



to decouple k and k' integrations in © and (ITU|) . after 
integration over k, k' we get 

(0 | H e i(r) | 0) = — — / d V . „ ■ ' (14) 



(r 2 + ry 2 ) e 

~ ~ AaTic f°° &r 2 n 2 

<0|«mW|0) = ~^-jf *7^+^)S- (15) 

Integration on 77 finally gives the well-known results (valid 
for r 7^ 0) 



(0 I H e z(r) I 0} 



23 a?ic 



(16) 



(4tt) 2 r 7 ' 

(0|^W|0) = -(4^^ (17) 

(22I [23l ] . The expressions above describe the local en- 
ergy densities of the electromagnetic field surrounding 
the point-like source. They are well defined everywhere, 
except at the origin r = 0, where they diverge. Before 
exploring in detail this singularity and show the existence 
of other singular terms at r = 0, it is of interest to dis- 
cuss some global (integral) property of the renormalized 
electromagnetic energy surrounding the dressed source. 
First, using ([2]) and (|6]), we immediately obtain 



(0 I H F I 0) = (0 I 



I 0) = 0. 



On the other hand, equations p^|) and (|17p imply 



(0 I {U el {v)+U m (r)) I 0) 



1 otht 



(18) 



(19) 



which diverges when integrated over all space because 
of the singular behavior of (fT6|) and (fl~7| at r = 0. 
This is at variance with the result (fT8|). because Hf = 
J d 3 r(H e i(r) + H m (r)). Thus there is a discrepancy be- 
tween the value of the total electromagnetic field energy 
calculated as the expectation value of the field Hamilto- 
nian Hp and as a space integral of the electromagnetic 
energy density (|TO1) . It appears that the average field en- 
ergy cannot be obtained from the spatial integration of 
(flU)) : also, the latter is divergent at the position of the 
atom. The mathematical origin of this difficulty is that 
the energy densities as given in (|16p and (|17[) have a sin- 
gularity at r = 0; we shall show in the next part of this 
Section that other singular terms are present and that 
their existence solves this (apparent) problem. This sin- 
gular behavior indeed prevents from exchanging the r and 
?7 integrations in the calculation of the total electric and 
magnetic energy densities [24| . To strengthen this obser- 
vation, we note that if we first perform integration over r 
and then that over r\ (that is in the opposite order of the 



r integral of (|19p ). the spatial integral of the total (elec- 
tric plus magnetic) energy density correctly vanishes, as 
expected from (TT5|) . In fact, from eqs. © and (fTTj) after 
integrations on k and k' , or from (TT4")) and (fT5j) . we have 



d 3 r(Q I H el (r) | 0) = f°° d V [°° dr 



3r 4 - 2rV + 3rj 4 hca f°° , 3tt 



(r 2 + r] 2 ) 6 



(20) 



d 3 r{0\H m (r)\0} = -^^ I drj I ,//-,-- 



8r 2 r] 



2, ,2 



Hca 



I" JO 

3tt 



dr] 



( r 2 + ^2)6 n 2 ^ ' 4r? 5 ' 

from which we immediately obtain 

d 3 r(0| (H e i(r)+H e i{r)) | 0) = 0, 



(21) 



(22) 



because the electric and magnetic contributions cancel 
each other, even if they are individually divergent. We 
can physically understand this result: the exponential 
function introduced in (fTBf plays the role of a cut-off 
function in the frequency integrations in (|9]) and (|10j) . and 
it removes the singular behavior of energy densities at r = 
before the integration over r\ is done. Exchanging the 
order of the integrations in 77 and r is somehow equivalent 
to remove the cut-off in the calculation, and this leads to 
the diverging result, essentially for the same reason that 
an ideal conducting boundary leads to diverging energy 
densities at the boundary- vacuum interface 14j . 

From the previous considerations, we are led to con- 
clude that the total field energy cannot be obtained by 
spatial integration of the energy density, without an ap- 
propriate prescription for dealing with the singularity at 
the origin. This observation suggests to investigate in 
more detail the singular behavior of the electromagnetic 
energy density. We will now show that if the singularity 
at r = is correctly evaluated and taken into account, 
the spatial integral of (0 | H e i(r) + H e i(r) \ 0) vanishes, 
consistently with (fTS)) . In order to do that, we first con- 
sider the electric and magnetic energy densities, following 
a procedure similar to that used in (24[. It involves the 
use of an exponential cut-off in the integrals © and (|10p , 
and letting the cut-off frequency to go to infinity after the 
frequency integrations (details are given in the Appendix 
A). After some algebra we get 

in i n, Mini ^ Ca [ 2 ^ 23 . 10 ... . 

~^'' (r) + 3^ y '' (r)+ li <5(ra)(r) }' (23) 

+ ^'' (r) -3^ ,5 ' , ' (r) -T^ (5M(r) }' (24) 
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where the superscript to the delta function indicates the 
order of its derivative with respect to r. The first term 
on the right-hand side of (f23| and (|24|) coincides with 
(HHJ) and (jTTJ) , respectively. The other terms, involving a 
delta function and its higher-order derivatives, take into 
account the singularity of the electric and magnetic en- 
ergy densities at r — 0. This singularity was not included 
in the calculation yielding (jTBf and (TT7)) . Our procedure 
outlined in the Appendix A of first introducing an expo- 
nential cut-off function, and taking the limit of the cut-off 
frequency to infinity after having done the frequency in- 
tegrals, has thus allowed us to obtain the correct and 
complete expression of the energy- density singularity at 
r = 0. 

We can now easily evaluate the total energy of the field. 
Summing up (|2"5|) and (IM1) , and integrating over all space 
using the propetries of the distributional derivatives of 
5— functions, we finally obtain 



d 3 r(0 | (H e i(r)+H 



mag 



(r))|0) = — J dr 



16 



16 
3^ 



48 
15H 



28 



15r 4 



S(r) \ = (25) 



(cancelation of terms diverging for r — > should be 
noted). As mentioned, this shows that, by a proper ac- 
count of the divergences at the atomic position, the to- 
tal field self-energy vanishes, consistently with the global 
self-energy as obtained in (fT8|) . This means that the to- 
tal (diverging) electromagnetic energy for r > is ex- 
actly canceled by the electromagnetic energy stored at 
r = 0, as clearly shown by eq. (|25p . This property is 
not obtained if the calculation of integrated energy den- 
sity is performed without a careful consideration of the 
singularity at r = 0, and shows that the renormalized 
field energy vanishes by virtue of a intriguing cancella- 
tion among diverging quantities. This fully restores con- 
sistency between the renormalized expectation value of 
Hp on the dressed ground-state and the integrated en- 
ergy density over all space. Although the total renormal- 
ized self-energy vanishes, its electric and magnetic com- 
ponents individually are however divergent at the origin 
for a point-like source. 

We wish to stress that a similar discrepancy has been 
discussed in the case of quantum fields with a bound- 
ary condition, for example for the electro mag netic field 
in presence of a perfectly conducting plate jlll.H4l| ; it has 
been shown that the existence of a singular surface energy 
density is necessary to remove the inconsistency between 
the total field energy and the integrated energy density 
[lij . We finally observe that the exponential cut-off we 
have introduced in the previous calculations is just a 
mathematical procedure that allows us to obtain the cor- 
rect form of the divergences of energy densities through 
a limit procedure. In the next Section we shall discuss 
how the assumption of a model of extended source may 
provide a natural physical cutoff, which removes all di- 
vergences in the energy densities of electromagnetic field. 



III. ENERGY DENSITIES AND 
FLUCTUATIONS OF THE ELECTROMAGNETIC 
FIELD NEAR AN EXTENDED SOURCE 

We now focus on the second main point of this paper, 
specifically the local and global electromagnetic energy 
densities in the case of a specific model of an extended 
source. We consider a neutral source of finite dimension, 
that we model as a collection of neutral point-like sources, 
smeared-out over a finite volume, with density p(t). This 
introduces a form-factor in the frequency integrations. 
Specific examples of form factors for transitions between 
levels of the hydrogen atom can be found in 26]. We 
assume the function p(r) normalized, so that 



d 3 iy>(r) = 1. 

A simple generalization of Hamiltonian (TT|) to the case of 
an isotropic source of finite dimensions is 

H = H A + ^- J d 3 r(d 2 (r)+b 2 (r)) 

kj Vj> v 
x J d 3 r'p(r')d k ^(r')) . (26) 

This Hamiltonian reduces to (fT]) for a point-like source 
when p(r) — > <5(r). In (|26| the interaction with the field of 
an elementary source at r' also depends on the presence 
of the elementary source at r". 

The density function p(r) appearing in the interac- 
tion Hamiltonian, acts as a regularizing factor of the 
source-field interaction: thus our extended-source model 
allows to remove all divergences in the renormalized val- 
ues of electric and magnetic energy densities, yielding 
also consistency between integrated averaged energy den- 
sities and average of the field Hamiltonian. Our model 
also permits to obtain the case of the point-like source 
discussed in the previous Section as a limit case. 

In order to evaluate the electromagnetic energy den- 
sity around the dressed source, we follow the same pro- 
cedure of the previous Section. We first calculate the 
dressed ground state of Hamiltonian ([26)) . Straightfor- 
ward perturbation theory up to first order in the atomic 
polarizability yields 



kj k'j' 



(fcfc') 1 / 2 . 
k + k' ' k '" k '' 



x II d 6 Td 6 r' p{r)p{r')e 



l\ — i(k-r+k'-r') 



<7,l kj -lk'i'>. (27) 



This expression can be used to evaluate the quantum 
averages (H e i{r)) and (% m (r)}. At lowest order in a{k) 
and subtracting the spatially homogeneous zero-point 
contributions (present also in absence of the source), we 
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get 



(Si*.(*>i5>....^) 2 E£ ( "' )I/a 

kj k'j' 



k + k' 



x(e kj -e k , r ) 2 (a(k) + a(k'))e^ +k '> R 

x / d 3 rp(r)e- ikr / dV ( o(r')e- 4k '- r ' + c.c. (28) 



for the renormalizcd electric energy density, and 



(0|^(R)|0) e .,=A(-) 



{kk') 



t\l/2 



k + k' 



kj k'j' 

x (em- • £k'j'){bkj ■ k'f)(a(k) + a(fc'))e l(k+k,) ' R 
x / d 3 r j o(r)e- jk r / d 3 r' j o(r')e- lk ' r ' + c.c. (29) 



for the renormalized magnetic energy density (subscript 
e.s. indicates extended source). Comparison with ([7]) and 
([5]) allows us to write these quantities in the more com- 
pact form 



(0 | H el (R) | 0) e . s . = i- / d 3 rp(r) / ,f'v' ,,( r' > 



((0 | d(R - r)d(R - r') | 0>„, 
1 



(30) 



(0 | H m (R) | 0) e . s . = — / d 3 rp(r) / dVp(r') 



((0 | b(R - r)b(R - r') | 0) p . s .) 



(31) 



(where the subscript p.s. indicates point-like source), 
which clearly shows the presence of interference between 
the contributions of the elementary sources which form 
the extended source. 

Assuming p(v) with a spherical symmetry, the form 
factor depends only on k, 



J d 3 rp(r)e- jk r = p(k). 



(32) 



Substituting (|32]l in (|30|) and (|3Tj) , in the continuum limit 
we obtain 

(0 I H el (R) | 0) e . s . = T7r~~3 / dk dk' (a(k) 

tD 71 " Jo Jo 

+ a(k')) {p(k)p{k')Q E {k, k\ R)-£—jj + c-c.J (33) 



and 



(0 | H m (R) | 0) e 



oc poo 



+ a(fc')) Uk)p(k')Q M (k, k', R)^^ + c.c.) .(34) 

These expressions show how the renormalized electric 
and magnetic energy densities explicitly depend on the 



structure of the source through its form factor p(k). In 
the limit p(r) — ► S(r), they reduce to the expressions for 
a point-like source obtained in the previous Section. 

The assumption of an extended source naturally in- 
troduces a physical cutoff (at a frequency related to the 
size of the source) , that eliminates the divergences in the 
field energy densities discussed in the previous Section. 
In fact, it is easy to see that the functions inside the 
integrals (|3"3")l and (jMJ) behave as a(k)p(k)k for large 
k. Thus choosing an appropriate form factor, decreas- 
ing as l/k a (a > 2), the integrals on the right-hand 
side of (|33[) and (|34|) converge. Furthermore, because 
the functions QE(M)(k, k', R) are continuous everywhere, 
both the electric and magnetic energy densities are well- 
defined everywhere. This also makes convergent their 
integral over all space. In other words, the singulari- 
ties present of the local energy densities in the case of 
point-like source are eliminated in our case of an extended 
source. 

Finally, we may calculate explicitly the spatial integral 
of the electromagnetic energy density. Now the integral 
over r can be safely exchanged with the integrals over k 
and k' due to the regularization introduced by the form 
factor; being 



d 3 r (Q E (k, k', r) - Q M (k, k' , r)) = 0, (35) 



we immediately get 

d 3 r(0 | (H el + U m ) | 0) e .». = , (36) 



as expected from eq. (|18|) . which is valid also for the 
extended source. This result confirms that because of 
the assumption of a source of finite dimension, there is 
not discrepancy between the total electromagnetic energy 
and the integrated local energy densities for any nonvan- 
ishing size of the source, making the procedure adopted 
fully equivalent to the introduction of a suitable cutoff 
function. 



IV. CONCLUSIONS 

In this paper we have analyzed the renormalized vac- 
uum electric and magnetic energy densities (or equiva- 
lently field fluctuations) near a field source, for example 
an atom or an electrically polarizable body, using an ef- 
fective Hamiltonian describing the atom-field interaction. 
The importance of considering these local energy densi- 
ties is related to many factors, for example they are re- 
lated to Casimir-Polder forces on polarizable bodies and 
they should also act as a source term for gravity. We have 
first concentrated our interest on the structure of the di- 
vergences of the energy densities at the source position 
(r = 0) in the case of a point-like source. We have found 
that the expressions of the electric and magnetic energy 
densities contain terms proportional to the Dirac delta 
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function and its derivatives evaluated at r = 0, which 
contain a finite amount of energy. We have shown that 
these singular terms in the energy densities are essential 
in order to have consistency between the renormalized 
vacuum expectation values of the field Hamiltonian and 
of the field energy density integrated over all space. We 
have also considered a model of an extended source of the 
electromagnetic field, which is smeared out over a finite 
volume of space, and shown that in this case the renor- 
malized field energy densities are finite at any point of 
space for any nonvanishing size of the source. Relation of 
our results with recent works about surface divergences 
of renormalized field energy density at the interface be- 
tween a conducting plate and the vacuum has been also 
discussed in detail. 
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Appendix A: Calculation of (H e i(r)} and (Hm(r)) 



drji / dkk 3 j 1 (kr)e- ( ~ f+v)k 



(A2) 



with 7 > and where we have used the relation (| 13)) . 
This integral can be easily evaluated; after some algebra 
we obtain 



J(r) = 64r 2 / drj 
32 



(7 + V) 2 



2 

r 7 



o V + ( 7 + ?7) 2 ] 6 

4 7 14 



7 



5(r 2 +7 2 ) 5 5 (r 2 +7 2 ) 4 15r 2 (r 2 + 7 2 ) 3 
77 77 



6r 4 (r 2 +7 2 ) 2 4 r 6 (r 2 + 7 2 ) 
7 1 T du 



4 r 6 ,/ 7 r 2 + /i 2 



(A3) 



where /i = 7 + 77. We now consider the limit 7 — > 0, 
equivalent to removing the exponential cut-off, after the 
frequency integrals. Using the Lorentzian representation 
of the Dirac delta function, it is easy to see that the 
terms appearing in the expression (| A3|) lead to the Dirac 
delta function and its derivatives when 7-^0. Thus, we 
obtain 



In this Appendix we outline the procedure that leads 
to the expressions ([23]) and (l24l) giving the energy-density 
singularity at r = 0. For simplicity, we first focus on the 
magnetic energy density (~H m (r)), given in (fT0|) . We need 
to calculate the following integral 



00 c oc 



k 3 k' 3 

I(r)= j dk J dk'nikr^k'r)^—^ 



(Al) 



In order to evaluate it, we introduce an exponential cut- 
off so that the integral becomes 



r°° k [i k Ki 
Hr) I dk dk'j 1 (kr)j 1 (k'r)f^ i e-^ k +^ 



(A4) 



where the superscript to the delta function indicates the 
order of its derivative with respect to r. Substituting in 
(fl"0|) and using (fl"2l) . we finally obtain expression (|24"|) . A 
similar procedure leads to expression (|23p for the electric 
energy density. 
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